Chiral effective field theory calculations of neutrino processes in dense matter 
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We calculate neutrino processes involving two nucleons at subnuclear densities using chiral effec- 
tive field theory. Shorter-range noncentral forces reduce the neutrino rates significantly compared 
with the one-pion exchange approximation currently used in supernova simulations. For densi- 
ties p < 10^*gcm~^, we find that neutrino rates are well constrained by nuclear interactions and 
nucleon-nucleon scattering data. As an application, we calculate the mean-square energy trans- 
fer in scattering of a neutrino from nucleons and find that collision processes and spin-dependent 
mean-field effects dominate over the energy transfer due to nucleon recoil. 
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Neutrino processes involving two nucleons (NN) play 
a special role in the physics of core-collapse supernovae 
and neutron stars: Neutrino-pair bremsstrahlung and ab- 
sorption, NN ^ NNvv, are key for the production of 
muon and tau neutrinos, for their spectra, and for equili- 
brating neutrino number densities [l|, 0, 0, 0; [Sj tSi] ■ Since 
neutrinos interact weakly, the rates for neutrino emis- 
sion, absorption and scattering are determined by the 
dynamic response functions of strongly-interacting mat- 
ter. Supernova explosions are most sensitive to neutrino 
processes near the protoneutron star, at subnuclear den- 
sities p < po/10 (po — 2.8 X lO^'^gcm"'^ being the sat- 
uration density), and where matter is neutron rich. In 
this regime, there exist to date no systematic calculations 
of neutrino rates that go beyond the one-pion exchange 
(OPE) approximation for the nucleon-nucleon interac- 
tion. In this paper, we present first results for neutrino 
processes in neutron matter based on chiral effective field 
theory (EFT). 

Noncentral contributions to strong interactions, due 
to tensor forces from pion exchanges and spin-orbit 
forces, are essential for the two-nucleon response. 
This follows from direct calculations of neutrino-pair 
bremsstrahlung 7] and from conservation laws |[8i] . In su- 
pernova and neutron star simulations, the standard rates 
for bremsstrahlungand absorption are based on the OPE 
approximation 0,l3|- This is a reasonable starting point, 
since it represents the long-range part of nuclear forces, 
and for neutron matter, it is the leading-order contribu- 
tion in chiral EFT However, for the relevant Fermi 
momenta kp ~ 1.0 fm~^ « 200 MeV, subleading noncen- 
tral contributions are crucial for reproducing NN scat- 
tering data 0. In this paper, we go beyond the OPE 
approximation and include contributions up to next-to- 
next- to-next-to- leading order (N'^LO) in chiral EFT. We 
find that shorter-range noncentral forces significantly re- 
duce the neutrino rates for all relevant densities. As an 
application, we calculate the mean-square energy trans- 
fer in scattering of a neutrino from nucleons and find that 



collision processes and spin-dependent mean-field effects 
are more important than nonzero momentum transfers to 
the nucleons. This establishes that the long- wavelength 
approximation used in the OPE calculations Q is rea- 
sonable. 

We follow the approach to neutrino processes in nu- 
cleon matter developed in Ref. [13], which is based 
on Landau's theory of Fermi liquids and consistently 
includes one-quasiparticle-quasihole pair states (corre- 
sponding to elastic scattering of neutrinos from nucle- 
ons) and two-quasiparticle-quasihole pair states, which 
are taken into account through the collision integral in 
the Landau transport equation for quasiparticles. Us- 
ing a relaxation time approximation, the transport equa- 
tion can be solved and this leads to a general form for 
the response functions [l3|- The spin response includes 
multiple-scattering effects, thereby taking into account 
the Landau-Pomeranchuk-Migdal (LPM) effect, and gen- 
eralizes earlier work to finite wavelengths. 
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FIG. 1: (Color online) Spin relaxation rate given by Co- of 
Eq. Q as a function of Fermi momentum fcp obtained from 
chiral EFT interactions of successively higher orders [3] ■ AH 
results are for m*/m = 1. 
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Basic formalism. In supernovae and neutron stars, the 
energy u and momenta q transferred by neutrinos to the 
system are small compared to the momenta of nucleons. 
In addition, we consider degenerate conditions, where 
the temperature is small compared to the Fermi energy, 
T/ep < 1/3. This is the regime in which Landau's the- 
ory of Fermi liquids is a reasonable first approximation, 
and includes the conditions under which the two-nucleon 
response is effective. We focus on axial current processes, 
since they dominate in many situations and are most af- 
fected by NN interactions. The corresponding neutrino 
rates are given by [if 



T{uj,q) = 27r n GpCl (3 -cos 61) 5a (t^,q) 



(1) 



multiplied by occupation probabilities for the initial neu- 
trino states and Pauli-blocking factors for final states. 
Here n denotes the neutron number density, Gp the Fermi 
coupling constant, Ca = —ga/"^ — —1.26/2 the axial- 
vector coupling for neutrons, and 9 is the angle between 
the initial and final neutrino momenta for scattering, or 
between the neutrino momentum and minus the antineu- 



trino momentum for neutrino-pair bremsstrahlung and 
absorption. The axial or spin dynamical structure factor 
is given by [11] 
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Imx^(tj,q) , 



(2) 



where x<y is the spin-density-spin-density response func- 
tion, and we use units with h — c = k^, = 1- The 
solution [l^ to the Landau transport equation leads to 
Xcr — N{0) Xa/{\ + GoXa), with density of states at the 
Fermi surface A^(0) = m*kp/Tr'^, nucleon effective mass 
TO*, and Landau parameter Gq for the spin-dependent 
part of the interaction. In the relaxation time approxi- 
mation, the response function (for Go — 0) is 



X„ = l- 



uj , f LU + i/Trj + vpq 
■ In' 



2vpq \LU + i/Ta- — v-pq 



(3) 



where denotes the spin relaxation time and vp — 
kp/m* is the Fermi velocity. For frequencies |a;| ^ l/rvr 
and q ^ kp, the spin relaxation rate is given by 



= G, + (c./2^)2] with G, = ^^^ ^ Tr ^,„^,(k,k')(Tf[(<Ti+^2f ,^.„..(-k,k')] 



fc=l,2,3 



, (4) 



where ^cti.<t2 (k, k') is the quasiparticle scattering ampli- 
tude multiplied by A^(0), k = pi— ps and k' = pi— p4 are 
the nucleon momentum transfers, and the average (. . .) 
is over the Fermi surface (for details see ^Oj). The com- 
mutator with the two-body spin operator demonstrates 
that only nonccntral interactions contribute. 



We treat the strong interaction in Born approximation 
and evaluate the spin trace for G^ in two-body spin space 
|s TOs) using a partial- wave expansion. For s = the spin 
trace vanishes, thus only s = 1 states and odd orbital 
angular momenta £,£',£,£' contribute, and we find 
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where j, j are total angular momenta, p is the magnitude 
of the relative momenta p = |pi — P2I/2 = |p3 — P4I/2, 
and (p|V^'^^|p) are partial- wave matrix elements of the 
strong interaction. We use a = \/2a + l and standard 
notation for Clebsch-Gordan, 3j and 6j symbols [13]. 

Results. In Fig. [1] we show the spin relaxation rate 



calculated from chiral EFT interactions up to N'^LO as 
a function of the Fermi momentum over a wide density 
range (fcp = 1.4 fm~^ corresponds to p = mfcp/(37r2) — 
1.6xl0^^gcm~'^) [i^]. The neutrino-pair bremsstrahlung 
and absorption rates are proportional to Ca to a good 
approximation when ^ I/t'o-, while for larger val- 
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FIG. 2: (Color online) The rms energy transfer AE in neu- 
trino scattering from neutrons as a function of initial neutrino 
energy E^, for two T-kp combinations. Results are shown for 
different Ca (OPE, N=^LO using Eq. (|9]), and a lower value), 
in the long- wavelength limit (5 = 0), all with Go ~ 0.8 [lq |. 
and in the absence of mean-field effects (Go — 0). 



ues of Co- they are suppressed by the LPJVI effect [l|, ■ 
The leading order (LO) contribution includes OPE as the 
only noncentral interaction, which provides the standard 
two-nucleon rates used in current supernova simulations. 
Our results based on chiral EFT interactions of succes- 
sively higher orders of Epelbaum et al. (EGM) [13] show 
that OPE significantly overestimates Co- for all relevant 
densities. The bands at next-to-leading order (NLO), 
N^LO, and N'^LO provide an estimate of the theoretical 
uncertainty, generated by varying the cutoff A as well as 
a spectral function cutoff in the irreducible 27r-exchange 
AsF Chiral EFT interactions at N^LO accurately 

reproduce low-energy NN scattering IJ, |l5| , and at this 
order, C^. is practically independent of the N'^LO po- 
tential for these densities. JVIost of the reduction of Co- 
occurs at the NLO level, which includes the leading 2tt- 
exchange tensor force (for Agp 00) 
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and shorter-range noncentral contact interactions 



(6) 



NLO 



C, nc 



Ci <Ti • k (T2 • k + C2 <Ti • k' cr2 • k' 

+ C3 I (<Ti - CT2) ■ (k X k') , (7) 



where w — ^/Am^^ + W and Ci are constrained by NN 
scattering data (S- and P- waves). ^2^^? added to OPE 
increases (to 0.29 MeV~^ at = 1.0 fm"^), but this 
is compensated by the repulsive shorter-range Vq^. 

Energy transfer. Figure [2] shows the impact of Co- on 
the root-mean-square (rms) energy transfer Ai? in scat- 
tering from nucleons. In the absence of Pauli blocking of 



> 

c 

O 



> 

e 

O 



0.25 

0.2 
0.15 

0.1 

0.05 


0.1 

0.05 - 







' I ' I ' I ' I ' I 
one-pion exchange 

n'lo 

V, , A=1.8-2.8 fm"' 

low k 

phase shifts 

m*/m=l 

H — I — I- 



low k 

fit Eq. (9) 




0.2 0.4 0.6 0, 



kp [fm" ] 



FIG. 3: (Color online) Combined results for the spin relax- 
ation rate given by Go-: bands for different chiral N^LO poten- 
tials [1^, tl5|] and chiral Viow fc interactions [ij and the result 
based on phase shifts. We also include the rate of Ref. 
for Viowfc with a density-dependent cutoff A = \f2kp. All 
results are for ra' jm = 1 (also for Ref. [13] )■ Lower panel: 
The band represents the Go- range in the upper panel based 
on nuclear interactions and NN phase shifts. The solid line 
includes second-order many-body and m* contributions, with 
m*/m « 1.12 — 0.17 (fcp/ fm~^) over this range p^ . 



final neutrino states, Ai? is given by 



J dpi [E, - Elf T{E, ^ El,p, - pI) 
JdplTiE,-El,p,-pl) 



(8) 



For the N'^LO results, the comparison with the long- 
wavelength approximation (g = 0) shows that the energy 
transfer due to nucleon recoil is smalL Collision pro- 
cesses are the major contributor to the energy transfer. 
This is also clear from the decreased AE with the lower 
Co- value. In addition, the comparison with the Gq = 
results demonstrates the increased energy transfer due to 
repulsive mean-field effects in the spin channel. 

Comparison of interactions. We have also carried out 
calculations for the Entem and Machleidt (EJVI) N^LO 
potentials with A = 500 and 600 MeV 15]. The EJVI 
500 MeV results overlap with the EGM N^LO band of 
Fig. [H The EM 600 MeV rate is larger due to the failure 
of the Born approximation in this case. We have also 
used the renormalization group (RG) [13] to evolve all 
N'^LO IJ, UM potentials to low-momentum interactions 
Viowk with A = 1.8-2.8 fm"^ (360-560 MeV) to extend 
the estimate of the theoretical uncertainty. The RG pre- 
serves the long-range pion exchanges and includes sub- 
leading contact interactions, so that NN scattering data 
are reproduced. The resulting Co- band is shown in the 
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upper panel of Fig. [3l This includes the rate based on 
the evolved EM 600 MeV potential. 

In the low-density limit, two-nucleon collisions domi- 
nate and NN phase shifts provide a model-independent 
result for the spin relaxation rate. In the near-degenerate 
case, the strong potential V in Eq. ^ is then replaced 
by the free-space K matrix, since Pauli blocking removes 
the imaginary parts from the loop integrals in the T- 
matrix equation. We use the standard expression of the 
K matrix in terms of empirical phase shifts and mixing 
angles j3| from the Nijmegen Partial Wave analysis (nn- 
online.org). In general, it is however unclear whether an 
expansion in terms of the free-space K matrix (or the T 
matrix) is reliable for kp > 0.1 fm~^ due to large scatter- 
ing lengths and Pauli blocking for near-degenerate con- 
ditions, unless the relevant parts of V are perturbative. 
Neutrino-pair bremsstrahlung neglecting mean-field and 
LPM effects has been calculated from the T matrix in 
Ref. For all densities in Fig. [31 we find the striking 
result that the chiral EFT and Mowfe rates obtained in 
Born approximation are close to those from phase shifts 
(with the K matrix). This demonstrates that, for these 
lower cutoffs, the noncentral part of the neutron-neutron 
amplitude is perturbative in the particle-particle channel. 

The upper panel of Fig. [3] combines our results based 
on nuclear interactions and NN phase shifts, and also in- 
cludes rates from Viowk interactions [l3| that extend to 
higher densities. At subnuclear densities p < 10^"* gcm~'^ 
(fcp < 1.2 fm"^), the spin response is well constrained and 
all results lie within a band, with a significantly reduced 
Co- compared to OPE. For kp > 1.5 fm~^, calculations of 
the equation of state show that low-momentum 3N inter- 
actions become important (lOl and should be included. 
In the lower panel of Fig. [3l we show a simple fit repre- 
senting our results. 
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(9) 



that, together with the spin response function given by 
Eqs. ^ and ([3]), can be used in astrophysical simulations. 

Many-body effects. Our results above do not include 
the effects of the nuclear medium on collisions. To give a 
sense of how these affect , we show in the lower panel 
of Fig. [3] rates that include second-order many-body con- 
tributions (particle-particle/hole-hole and particle- hole) 
and self-energy effects through the effective mass m* (ac- 



cording to Eq. (O, Ca 



') [10]. Both particlc-hole 



and m* effects reduce the spin relaxation rate. The for- 
mer is driven by second-order particle-hole mixing of ten- 
sor with strong central interactions 2l|. This demon- 
strates the need to study in greater detail the influence 
of many-body effects on collisions (see also Ref. ^22j ) . 

In summary, we have presented the first calculations 
of neutrino processes in supernovae based on chiral EFT. 
Our N^LO results over the important density range p < 



10"'^^gcm~^ represent a significant advance beyond the 
OPE rates currently used in simulations. Shorter-range 
noncentral forces reduce the rates significantly. For densi- 
ties p < lO^'^ gcm~^, the spin response is well constrained 
by nuclear interactions and NN scattering data. Future 
work will include neutron-proton mixtures and charged 
currents, to systematically improve the neutrino physics 
input for astrophysics. 
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